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Abstract 

We show that the standard Calogero Lax matrix can be interpreted as a function on 
the fuzzy sphere and the Avan-Lalon r-matrix as a function on the direct product of two 
fuzzy spheres. We calculate the limiting Lax function and r-function when the fuzzy sphere 
tends to the ordinary sphere and we show that they define an integrable model interpreted 
as a large N Calogero model by Bordemann, Hoppe and Theisen. 
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1 Introduction 



The formalism based on the concepts of the Lax matrix and the Lax pair turns out to be very 
efficient in a study of a large class of classical completely integrable models [U [7J [SJ |33]. Briefly 
speaking, the Lax pair is a pair of matrix valued functions L, M on the phase space P of the 
integrable model such that the equations of motion of the latter can be written in terms of the 
matrix commutator 

(1-1) 

Moreover, the Lax matrix L is required to generate integrals of motion in Poisson involution, 
i.e. it should fulfil a condition 

{tr L m , tr L n } P = 0, Vm,n. (1.2) 

Bordemann, Hoppe and Theisen [10] had an idea to generalize this formalism by replacing the 
matrices L, M by functions on some two-dimensional auxiliary compact symplectic manifold B 
and the commutator by the auxiliary Poisson bracket {., .}#: 

j t L = {L ) M} B . (1.3) 
Moreover, they imposed the following condition on L 



uL m , uL n \ = 0, Vm,n, (1.4) 

B J B ) p 

where u denotes the auxiliary symplectic form on B. 

Starting with a suitable ansatz for the Lax function, Bordemann, Hoppe and Theisen [TU] have 
identified several field theories in 1+1 dimensions which are integrable via their novel mechanism. 
Among them they singled out two models the classical actions of which read respectively 

S t = \ j dt j rf(7^(a) 2 -a 2 expg'(a)^, (1.5) 

Sc = iJ dt J\^ ? -j^y (1 . 6) 

Here the dependence of the dynamical field q(a) on the time t is tacitly understood, a is a 
coupling constant, the "dot" and "prime" denote the time and the space derivatives, respectively, 
and boundary conditions on q(a) will be specified in Section 3. 

It is instructive to compare the field theories fll.5p . fll.6p with the integrable iV-particle Toda 
and Calogero systems: 

StN = ^ / dt ( X] < li~ K N ex P(^-9i+i)), (1-7) 

J M<i<AT l<i<N ' 

J K l<i<N l<i^j<N w ^ v 7 
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The similarity between field theoretical and iV-particle expressions was not the only reason 
why Bordemann, Hoppe and Theisen interpreted the models fll.5p . fll.6p as "infinite dimensional 
analogues of iV-particle Toda- and Calogero systems" . In fact, although they did not derive these 
field theories as large N limits of the iV-particle models, they argued that it should be possible 
to do it because the auxiliary Poisson algebra of functions on the two dimensional compact 
symplectic manifold B can be viewed as infinite-dimensional version gl(oo) of the algebra gl(N) 
of finite size N x N- matrices [27]. The parameter N/2 should be then interpreted as the Planck 
constant of an auxiliary quantization converting the auxiliary Poisson algebra into the matrix 
algebra and the integrable field theories fll.5p . fll. 61) could be therefore called "doubly classical" 
or "totally classical". 

The suggestion of [10] to obtain the Lax functions of the totally classical field models fll.5p . fll.6p 
as a large iV limit of the known Lax matrices of the iV-particle models fll.7p .f jL8|) was then 
realized by Hoppe in the Toda case [29]. However, to our best knowledge, a similar result was 
not obtained for the Calogero model. We wish to stress that the similarity between the field 
theoretical model and the finite iV-particle theory is more evident for the Toda model than for 
the Calogero one. Indeed, the Toda field theory (11. 5p is local in the sense that its Lagrangian 
contains only the first derivative q'{cr) and the discrete theory (11. 7p is also local in the sense that 
the potential contains only the interaction of the nearest neighbours. In reality, the iV-particle 
Toda model (11.71) is nothing but a straightforward discretization of the Toda field theory (11.51) . 

In the Calogero case the situation is much more subtle, because the field theory (11.61) is local 
but the iV-particle Calogero model (11 .8p is non-local. Indeed, the Calogero model (11 .8p is not 
a discretization of the field theory (11. 6p . On the top of it, there exists a non-local field theory 
referred to as the large N Calogero model in the literature [U UT] the discretization of which 
does give the iV-particle Calogero model (11.81) . Would the existence of two different large N 
limits of the same finite N theory indicate a contradiction? A priori no, since much depends 
in which context the large N limit is taken, what other parameters get fixed while taking the 
limit etc. [51 [HI I2H [23 ED E21 ESI EZ] • We shall actually justify the interpretation of the model 
(11.61) as one particular large N limit by appropriately tuning the coupling constant and the 
Darboux symplectic Qn form of the iV-particle Calogero model (11. 8p . Speaking more precisely, 
we shall prove that the discrete non-locality gets eliminated and the finite N Calogero model 
approaches the totally classical field theory (11. 6p if we make both and Qn proportional to 
1/N. This will constitute the first result of the present paper. 

In the original article [10], the fundamental involutivity relation (11. 4p was obtained by a direct 
computation starting from an appropriate ansatz for the Lax function. In the finite N case, the 
involutivity (11. 2p can be proved most easily by using the concept of the classical r-matrix. We 
recall that the classical r-matrix [7J, |8] is a map r(N) u from the phase space of the model into 
the direct product of two copies of the matrix algebra such that a fundamental compatibility 
relation with the Lax matrix L(N) holds: 

{L(N) U L(N) 2 } P = -^[r(A0i 2 ,Li(A0] + ^ [r(N) 21 , L 2 (N)), (1.9) 
where L(N)i = L(N) <S> 1 and L(N) 2 = 1 <S> L(N). It is then very easy to show that (I1.9P implies 

In the totally classical context, Hoppe [28J has suggested to introduce a concept of an r-function 
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on P x B x B, which would verify an analogue of (II. 9ft : 



{L{z), L(w)} P = {r{z, w), L(z)} B - {r(w, z), L(w)} B . (1.10) 

Here L is the Lax function, the dependence of L and r on the coordinates of P is tacitly 
understood in (ll.lOp and the letter z or w stands for some parametrization of the auxiliary 
symplectic manifold B. As in the finite N case it is easy to show that the involutivity relation 
( II. 4p is the consequence of (11.101) . 

In the case of the totally classical Calogero model, Hoppe started to look for the r-function by 
choosing an appropriate ansatz which he then substituted into the condition (11.101) . Remarkably, 
the resulting object which he has found is not quite an r-function but rather an r-distribution 
in the variables parametrizing B x B\ This may look surprising at a first sight: how the large 
iV-limit may convert matrices into distributions? As we shall see, this is indeed the case and the 
r-distribution does arise in the large N limit. The derivation of Hoppe's r-distribution directly 
from the A van- Talon r-matrix [6] of the iV-particle Calogero model constitutes the second result 
of the present article. 

In this paper we choose the ordinary sphere as the auxiliary symplectic manifold B and in 
Section 2 we review its quantization called the fuzzy sphere [27, 35J. Then in Sections 3 and 4 
we prove that the Lax matrix and r-matrix of the iV-particle Calogero model are respectively 
the fuzzy quantizations of the Lax function [10J and of the r-distribution of the totally 
classical Calogero model (11.61) . We finish by a short outlook. 



2 The fuzzy sphere 



The ordinary symplectic sphere S 2 is a surface embedded in IR 3 defined in Carthesian coordinates 

xi, 22,^3 as 

x\ + xl + xl = l. (2.11) 

We parametrize it as 



X\ = Vl — <J 2 cos 0, X2 = Vl — cr 2 sin 0, x% = cr, a G [—1, 1], G [— 7r, 7r]. (2.12) 
The standard round symplectic form on S 2 is then given by 



u = --e jk iXjdx k A dxi 



= daAd(j). (2.13) 

s 2 



The result of the quantization of the symplectic manifold (S 2 , u) is known as "the fuzzy sphere" 
[2Z1E5]- The linear SO (3)-equivariant quantization map Qn associates to smooth functions / 
on S* 2 sequences of iV x iV-matrices Q,N{f) which are called the quantized or fuzzy functions. 
We shall not need an explicit formula for the quantization map Qn but we do need three basic 
properties of it: 

" 2 



QN(f)Q N (g) = QN[fg + 0[-)), (2.14) 



[Q N (f), Qn(9)} = Qn [ ijj{f,g} B + O ( — ) ) , (2.15) 
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Obviously the parameter 2/N plays the role of the auxiliary Planck constant. 

To give a flavor, what the map Qn is about, let us make explicit the quantized versions of the 
functions 23, X\ ± ix2 defined in (j2.12|) : 



Q N (xz)ij = -y===(N + 1 - 2j)S ij} Q N {x l + ix 2 ) t] = ^—— ^{j - 1)(N — j + l)^_i 

(2.17) 

and Qn(xi — 1^2) is the Hermitian-conjugated matrix Qn(xi + ix 2 Y . In particular, it is then 
easy to verify that it holds the emblematic fuzzy sphere relation 

Qn{xi) 2 + Qn(x 2 ) 2 + Q N (x 3 ) 2 = 1 N , (2.18) 

where ljv stands for the unit A" x A^-matrix. 

So far we have learned that every smooth function on the sphere gives rise to a sequence of 
iV x iV-matrices. It is perhaps less known that appropriate sequences of iV x iV-matrices may 
represent quantizations of not just smooth functions, but also of singular functions and/or 
distributions on S 2 . As an example particularly relevant for the present paper, we now describe 
the fuzzy vortices. For that, consider the sequence of N x N matrices V(N) defined by their 
matrix elements: 

HA'),, ,. z,j = l,...,iV. (2.19) 

In words: V(AQ is the Jordan block with zeros on the principal diagonal. We immediately check 
that 



Q N {x x + ix 2 ) = \jl- ^=±Q N ( X3 )^jl + ^=±Q N ( X3 ) V (N) (2.20) 

and we note that the diagonal matrices of which we take the square roots have all eigenvalues 
strictly positive. In the limit A^ — > 00, we then obtain 



X! + ix 2 = yjl-xi Kc, (2.21) 

which together with f)2.12p leads to an identification of Voo with the ordinary sphere vortex 
configuration e 1 ^ and to an extension of the quantization map to this vortex configuration by 
setting 

Q N (j*) = V{N). (2.22) 

Another relevant quantized singular object on the sphere is characterized by the following se- 
quence of N x N matrices K(N): 

N 

K(N) :=J2 E kh (2.23) 

k,l 

where Em is the elementary matrix with 1 in k th row and I th column and everywhere else. 
Note that K(N) are matrices with all elements equal to 1. For a finite N, we observe 

K(N) = (V(N)Y' 1 + - + (V(N)^ 2 + V(NY + 1 n + V(N) + V(N) 2 + ... + ViNf' 1 (2.24) 
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which gives from f |2.22j) 



K 00 = \im N ^ 00 Q N ( V = 2tt5(0). (2.25) 



Said differently, we extend the quantization map to the delta function by setting 

2irQ N {8{<f>)) = K(N). (2.26) 

Now we use the same symbol Qn for the quantization of the direct product S 2 x S 2 . We can 
argue that 

/ \ N N 

Q N I - tr 2 )5(0i - 2 ) J = j- J2 E ki ® (2-27) 

^ ' 71 k,i 

Indeed, the 5-function 8{a\ — <j 2 )5(4>i — 2 ) on S 12 x S 2 is characterized by the property 



i 



da 2 / tZ0 2 5(<Ti - cr 2 )5(0! - 4> 2 )f(cr 2 , 2 ) = / o; 2 5(cri - cx 2 )5(0i - (f> 2 )f(cr 2 , 2 ) = /(ci, 0i), 



(2.28) 

where /(cr 2 , 2 ) is an arbitrary smooth function on S 2 . Following f)2.16p . the quantized version 
of is 

^tr 2 ( Q,v(<f(<Ti - <r 2 )<K0i - 2 ))(1 JV <8) Qiv(/))1 = Qiv(/)- (2.29) 



N 

On the other hand, for the N x iV-matrix Qn(J) it obviously holds 

tr 2 r(^^®E Zfc )(l iV g)Q iV (/)) > ) =Q N (f). (2.30) 
Comparing the last two equalities leads to the identification (I2.27p . 

Finally, we shall need also the fuzzy version of the "diagonal" delta function 5(a\ — a 2 ) on 
S 2 x S 2 . Since the fuzzy version of a ^-independent function on S 2 is a diagonal matrix, the 
fuzzy version of 8{a\ — a 2 ) must be given as a sum of direct products of diagonal matrices. To 
find out which diagonal matrices appear in those direct products we look for a fuzzy analogue 
of the following identity: 

J da 2 J d(j} 2 S(a 1 - a 2 )f{a 2 ) = J u 2 5(a t - a 2 )f{a 2 ) = /(o"i). (2.31) 

It obviously reads 

| ? tr 2 (g iV (% 1 - ( 7 2 ))(l JV (8)g J v(/)) > l =QnU)- (2-32) 



N 

We have also the following matrix identity: 

N 



ti 2 (^rE kk ®E kk )(l N ®Q N (f))J =Q N {f), (2.33) 



which holds for the diagonal matrix Qn{J')- Comparing the last two equalities leads to the 
desired quantization formula 



N N 



1 V v — "\ 

Qiv(5(ai - a 2 )) = — E kk ® E kk . (2.34) 



3 Large N limit of the Calogero Lax matrix 

Recall first the standard formula for the Lax matrix L(N) of the Calogero model 

L(N) l3 = pAj + (1 - 1 < i,j, < N, (3.35) 

Qi Qj 

where pj, qj are the Darboux coordinates on the phase space of the model. In what follows, we 
normalize the Darboux symplectic structure Qn and the coupling constant kn as 

2 c 
fl N = jrdqj A dpj, k n = — . (3.36) 

We now introduce important diagonal N x iV-matrices R{N) and P(N) as follows 

R(N)ij := q^, P{N) ij := p^, i,j = l,...,N (3.37) 

and we naturally interpret R(N) as QAr(g(cr)) and P(N) as Qn(p(&))- Speaking more precisely, 
we have from (12.171) 

Qj = QN{q{a))jj = q { ^JjfT^ =f) ( 3 ' 38 ) 

and similarly for pj. 

Few comments are in order to justify the interpretation of R(N) as Qn(q(&)) and P{N) as 
Qn(p(ct))- First of all, q(cr),p(cr) are viewed as the functional coordinates of the field theoretical 
phase space of the totally classical Calogero model (11. 6p with the Darboux Poisson bracket 

{p(o-i),q(cr 2 )}p = 5^ - a 2 ), (3.39) 

but at the same time they are viewed as ^-independent functions on the sphere. We note also 
that for j running through the set l,...,iV the argument ^7^==^ in (I3.38P runs equidistantly 
through the interval [—1,1] which is indeed the domain of definition of the function q(o~). More- 
over, for neighbouring j and j + 1 the distance of the arguments is -r=g 1 which for large iV 
is nothing but the S-Planck constant. We thus observe that the phase space P of the totally 
classical Calogero model (II .6p gets drastically shrunk by the fuzzification. Indeed, among all 
functions q(o~),p(o~) which constitute points in P, only functions constant on the equidistant 
intervals of the lenghts survrve the fuzzification and form a 2iV-dimensional phase space 

P{N) of the iV-particle Calogero model. Moreover, we can derive from the Darboux Poisson 
bracket (I3.39P on the phase space P the Poisson bracket on the 2iV-dimensional phase space 
P{N) parametrized by qj,Pj, j = 1,...,N. To do that, we consider a function T(o~) on the 
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sphere and its diagonal quantization Q N {T). From the Darboux Poisson structure (13.391) . we 
obtain 

f uTqj =T(a). (3.40) 
Now from the property (12.161) we see that the fuzzification of the formula (I3.40p must give 



Pi, ^tr (Q N (T)R(N)) \ = Q N (T) U . (3.41) 

iV 1 P(N) 



We thus infer 

N 

{Pi,qj}p(N) = -j s v ( 3 - 42 ) 

which is consistent with the first formula in (I3.36|) . 

Let us verify that the iV-particle Calogero Lax matrix L(N) f!3.35j) has as a large iV-limit some 
classical Lax observable L. Said differently, we shall show that L(N) can be interpreted as the 
fuzzification Qn(L) of some function L on the ordinary sphere. 

It is now easy to verify with the help of (12.261) that 

[R(N),L(N)} = w N (K(N) - 1 N ) = ^Q N {2^5{(j>) - 1), (3.43) 

where the matrix K(N) was defined in (12.231) . If it is true that L(N) can be identified with 
Qn(L) for some L then we could rewrite (I3.43P as 

[Q N (q(a)), Q N (L)} = -^(2^(0) - 1). (3.44) 

which, due to the fundamental quantization property (12. 15 j) . would lead to 

{q(a), L(a, 0)} B = -q'(a)d^L = |(2tt5(0) - 1). (3.45) 

And indeed! The differential condition (I3.45P has the following obvious solution 

L(a, <P) = p(a) + — ^(0 - vrsign(0)). (3.46) 
2p'(<7) 

which coincide^] for c = 2y/3a/n with the Lax function of the model (ll.6p as found in [10] . 
Moreover, the condition (I3.45P determines the function L almost unambiguously. The only 
ambiguity consists in adding a 0-independent function to the solution (I3.46p . however, this 
ambiguity is fixed by the fact that the diagonal term in the Lax matrix L(N) must be 
equal to Qat(p(o")). Thus we have justified the interpretation [TU] of the totally classical Calogero 
model as the large iV-limit of the iV-particle Calogero model. 

Remark: It is important to stress that, as it stands, the formula (I3.46P defines only a func- 
tion on a subset of S* 2 covered by the coordinate chart (a, (f>) and we need the classical Lax 



1 The comparison of (|3.46[) with the result of [TU] must take into account the range of the parameter <f> which 
is [— 7r, 7r] in our paper and it was [0, 2ir\ in [TUJ. We have opted for the different range in order to stress that the 
Lax function L is discontinuous. This fact is indeed less visible in |10) since the discontinuity occurs precisely 
on the boundaries of the range [0, 2n]. 
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observable everywhere on S 2 . Fortunately, the function — 7rsign(0) smoothly extends to the 
anti-Greenwich meridian <fi = ±n and we can also extend L^^p, p; a, 0) to the poles a = ±1 if 
we impose the following boundary conditions 

lmv^ ± i// (cr) = +oo. (3.47) 

It can be also checked that these somewhat exotic boundary conditions are consistent with the 
dynamics of the model fll.fip since they make sure that there is no flow of energy through the 
boundaries a = ±1. 



4 Large N limit of the A van- Talon r- matrix 



Consider an r-matrix r(N) given by the formula 



N . 1 N i 

r{N) 12 = V — - — E kl ®E lk + -\2 — - — E kk ® {E kl - E lk ) 



kj4 



(4.48) 



Using the Poisson brackets (13.421) . it is a matter of straightforward computation to verify that 
the r-matrix (14.48)) and the Lax matrix L(N) satisfy the following crucial property (cf. fll.9D ) 



{L(N) u L(N) 2 } P{ir) 



~[r(N)i2,Li(N)} + 1 -^-[r(N)2i,L 2 (N)}. 



(4.49) 



We recall that it is this property ( 14.49)) that guarantees the P(iV)-Poisson commutativity of 
the traces tr L(N) n . For a computational convenience, the matrix r(N) slightly differs from the 
r-matrix proposed by Avan and Talon in [6]; in fact, r(N) is just another element of a moduli 
space [20] of all r-matrices verifying (14.491) . 

Recalling the matrix R(N) introduced in the previous section it is easy to check that 



[R(N) u r(N) 12 ] 



N 



E, 



Ik; 



[R(N) 2 ,r(N) 



12 



N 
k,l 



E 



ik 



N 



mm ® I'-min ) • I 1 



N 



N 



(4.50) 



(4.51) 



where as usual R(N)i = R(N) £g> 1, R(N) 2 = 1 <8)R(N) and [., .]+ stands for an anticommutator. 
Now writing R(N) as Qn(q(ct)), using the formulae (I2.26IUI2.27I) and (12.341) and supposing that 
r(N) can be written as a fuzzification Qn{j") of some function r on S 2 x 5* 2 , the equations (14.50)) 
and (I4.5ip can be rewritten as 



[Q N (q^))i,Q N (r)} 



2i 



Q N 2'k5(o 1 - a 2 )5((j)i - 4> 2 ) - 5(<T! - (T 2 , 



(4.52) 



2i 



[Qjv(?(o-))2, QN(r)] = JjQn [ 2n5(a 1 - <r 2 )5(0i - < > 2 



l 

N 



QN{S{(Tl-(T2)),QN{2'K6{ ( f> 2 )) 

(4.53) 
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Using the fundamental quantization properties f !2.14j) and f !2. 15j) . we deduce from (I4.52p and 
(EgD that r must fulfil 



{g(cri),r(<Ti,0i;cr2,02)}B 
{9(^2), »"(o r i,0i; o- 2 , 2 )}b = 



-g'(<n)<V = <J(<r a - a 2 )(l - 2tt5(0i - <j> 2 )). 
q'i^dfcr = 2n5(a 1 - a 2 )(5((f> 1 - <f> 2 ) ~ 5 (fa)). 



(4.54) 
(4.55) 



The relations (I4.54p and f 14 . 5 5 p are indeed verified by the following distribution 

r(oi, fa; (7 2 , 2 ) = —^8{a x - a 2 ) [E(fa - fa) + E(fa)) . (4.56) 

Here E is viewed as a 27r-periodic function on the whole real axis M. which is given by the 
expression — 7rsign(0) when restricted to the interval [— 7r, 7r]: 



The distribution (I4.56P does coincide with the one found in [28], moreover, the first order 
differential conditions (I4.54p and (I4.55P determine r almost unambiguously. The only ambiguity 
consists in adding a fa, ^-independent function to the solution (14.561) of (14.541) and (14.551) . 
however, such function would be a large N limit of a bi-diagonal term which is absent in r(N). 
We thus conclude that the classical Yang-Baxter observable r (14.561) of the totally classical 
Calogero model (II. 6p is indeed the large N- limit of the A van- Talon r-matrix r(iV)i 2 . 

5 Outlook 

As it is well-known, the Calogero model has attracted a lot of attention in pure mathematics 
and in mathematical physics (for reviews see [3l[T8|,[T9l[39lll0lll3llllj) since it appears in a large 



variety of contexts El HH 021 031 US1 HZl 1221 1231 1313 IBSl B21 SSI SSI SZl SSI SSI 



ranging from condensed matter physics, higher spin algebras, two-dimensional QCD or fluid 
dynamics to microscopic description of black holes etc. We expect a relevance of the totally 
classical Calogero model in many of these contexts. Speaking more generally, it would be 
interesting to construct the totally classical models corresponding to various trigonometric, 
elliptic or even relativistic deformations of the Calogero one. This would certainly contribute 
to a better understanding of the duality properties of the integrable systems. One can also look 
for gauge transformations of the r-functions in order to get rid of the q dependence of r (14.561) . 
much in the spirit of [20] which does that for the standard r-matrices. 
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